Discrete-time quadrature feedback cooling of a radio-frequency mechanical resonator 
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We have employed a feedback cooling scheme, which combines high-frequency mixing with digital 
signal processing. The frequency and damping rate of a 2 MHz micromechanical resonator embedded 
in a dc SQUID are adjusted with the feedback, and active cooling to a temperature of 14.3 mK 
is demonstrated. This technique can be applied to GHz resonators and allows for flexible control 
strategies. 



Mechanical systems in the quantum regime |l] [| can 
be used to answer fundamental questions about quantum 
measurement, decoherence, and the validity of quantum 
mechanics in macroscopic objects. This requires a me- 
chanical resonator which is cooled to such a low tem- 
perature that it is in its ground state for most of the 
time. In the past few years tremendous progress has 
been made in actively cooling resonators p), mainly by 
using sideband cooling |[ |J and active feedback cool- 
ing ||~||. The latter technique has mainly been applied 
to low-frequency (kHz) resonators combined with opti- 
cal detection. The largest cooling factors have been ob- 
tained using velocity-proportional feedback, i.e., by feed- 
ing back the differentiated displacement signal. However, 
at higher frequencies, delays in the feedback system seri- 
ously degrade the cooling performance. Here, we demon- 
strate a feedback cooling technique || with a nearly un- 
limited bandwidth, based on fast digital signal processing 
(DSP) in combination with single-sideband mixing. A 2 
MHz micromechanical resonator with inductive readout 
is cooled to 14.3 mK using this scheme. 

Figure |l|a shows the device, which consists of a dc 
SQUID with a part of its loop suspended. This forms a 
50 /im long flcxural resonator with its fundamental mode 
around /o ~ 2 MHz. The chip is glued onto a piezo el- 
ement for feedback and actuation, and cooled in a di- 
lution refrigerator with a minimum bath temperature of 
15 mK. By applying an in-plane magnetic field B (green), 
a displacement of the beam u changes the amount of flux 
through the dc SQUID loop. When a bias current Ib 
is applied, a change in flux results in a change in the 
SQUID voltage V. This way, the dc SQUID is a sensi- 
tive displacement detector Jl0| . In all measurements pre- 
sented here, the same working point for the dc SQUID is 
used, to avoid backaction-induced changes in frequency 
and damping pd[ |. 

The thermal noise of the resonator is used to calibrate 
the dc SQUID detector. Figure [|b shows the displace- 
ment noise spectrum S uu measured at two different cryo- 
stat temperatures T. The thermal motion of the res- 
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FIG. 1: (a) Schematic overview of the SQUID detector (red) 
with the integrated flexural resonator, (b) Displacement noise 
spectra without feedback, (c) The resonator temperature ex- 
tracted from the thermal noise spectra plotted against the 
mixing chamber temperature, (d) Generic linear system rep- 
resentation jl2| of feedback cooling, (e) The feedback filter 
consists of a digital signal processor with a single-sideband 
mixer at the input and output. 



onator shows up as a peak on top of the imprecision 
noise floor S UnUn . The cryogenic-amplifier-limited dis- 
placement noise is 2fm/\/Hz at T = 15 mK. The area 
under the peak is the amplitude of the Brownian mo- 
tion of the resonator squared. When the temperature of 
the refrigerator is increased to T = 0.3 K, the spectrum 
changes: Firstly, the noise floor is higher due to a de- 
crease in the transduction, dV/du, as the critical current 
decreases with increasing temperature [|l0| . Secondly, the 
peak is higher and wider (the intrinsic damping rate 70 
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increases with temperature) , indicating that the thermal 
motion is larger at higher temperatures. The resonator 
temperature T = k (u 2 )/kB (ko = 110 N/m is the spring 
constant) is plotted against T in Fig. 0c: It follows the 
cryostat temperature for T > 50 mK (solid line) and sat- 
urates below this value. 

To further lower the resonator temperature, active 
feedback is employed, where the displacement of the res- 
onator is fed back to it to damp its thermal motion. Fig. 
|l|d illustrates the generic process ||: The thermal force 
noise F t h = fco/t/i drives the resonator whose response is 
Hr = Z 7(/o - P+ */7o/27r). Note, that f th and the 
other signals are scaled to have the unit of position. The 
force results in a displacement which is measured by the 
dc SQUID detector, and imprecision noise u n is added to 
its output v. This signal is fed to the feedback filter with 
transfer function Gfb- The actuation a is multiplied by 
A, which consists of the SQUID transduction, an atten- 
uation (-40 dB) and the piezo responsivity. Finally, the 
resulting piezo displacement u p exerts an inertial force on 
the resonator. Note, that in practise crosstalk (X) exists 
between the applied feedback and the detector output, 
which modifies the system response. 

To fully characterize the linear system an ac signal is 
applied to a (see Fig. |l|d) and the response at v is mea- 
sured at the same frequency, while sweeping the driving 
frequency across the resonance. In this case, the feed- 
back Gfb is disabled. From this network-analyzer mea- 
surements the elements A = 1.94 • 10~ 4 exp(— 0.73i), and 
X = 0.26 exp(2.56z) of the linear systems are obtained as 
well as the parameters of Hr: fo and 70. The non-zero 
phase of A is due to the time it takes for the signal to 
travel through the whole system. If an analog differentia- 
tor would be used for Gfb, this delay causes the feedback 
to not be purely velocity proportional thus degrading the 
cooling performance. The DSP-based feedback presented 
here can compensate for this effect as demonstrated be- 
low. 

Our implementation of the feedback filter Gfb is shown 
in Fig. |l|e. The high-frequency input signal v is split 
and both branches are mixed with local oscillator (LO) 
signals with a 90° phase difference between them. This 
IQ mixer gives both quadratures v s and v c of the input 
signal. The LO frequency is f LO = 2.0492 MHz so the 
down- mixed signals oscillate at Jr—Jlo = 8.9 kHz. They 
are digitized and the DSP (Adwin Pro II at a sampling 
rate f s — 820 kS / s) applies the following transformation 
to the input signals to generate two output signals a c and 
a s : 

( a c \ _ f cosOfb - sm6 fb \ f v c \ . . 
\a s J g f b { sm8f b cosfl/i, J { v s J [L) 

These quadratures are then up-converted by the LO fre- 
quency with a second IQ mixer. The final result is a 
signal a at the original frequency that is phase-shifted by 
the feedback phase 6fb and multiplied by the feedback 
gain gfb, i.e. Gfb = g fb exp(i# fb) ■ The only frequency 
requirement for this mixing scheme is that the quadra- 



tures do not change faster than the sampling rate, which 
is equivalent to 7^/27r < is/2. The operation is thus not 
limited to resonators with frequencies within the band- 
width of the DSP, allowing feedback cooling of radio and 
microwave frequency resonators. Note, that the impreci- 
sion noise floor is still determined by the cryogenic ampli- 
fier; the contributions from the mixers and discretization 
are negligible. 
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FIG. 2: Feedback phase dependence of the resonator fre- 
quency (a), damping rate (b), and resonator temperature (c) 
for gft = 0.1. The dashed lines indicate their measured val- 
ues in the absence of feedback, (fo, 70, and To resp.); the 
solid line is the phase dependence calculated using indepen- 
dent measurements. 

The feedback modifies the resonator response from Hr 
to its closed-loop form H' R Q: 

P 

H ' R = p-P + iflo/2ir-pG> fb A> (2) 

where G' fb = g' fb exp(i9' fb ) = G fb /{1 - XG fb ) is the 
feedback filter modified by the crosstalk. The real part 
of G'j b A modifies the resonance frequency from /o to 
fa ~ /o(l — R- e [G^- b A]/2), whereas the imaginary part 
changes the damping from 70 to 7^ 70 — 2tt f lm[G'j: b A]. 
Both the frequency shift and the change in damping de- 
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pend periodically on the phase of G'f b A and the maxi- 
mum frequency shift is half the maximum damping rate 
change. 

In order to achieve optimal cooling the feedback phase 
is varied for a fixed feedback gain as shown in Fig. [|. 
The feedback gain is chosen sufficiently small so that 
G'j b ~ Gfb- At every point the thermal noise spectra are 
measured and fitted to obtain the resonance frequency 
(top), the damping rate (center), and the resonator tem- 
perature (bottom). The resonance frequency and damp- 
ing rate show the expected sinusoidal dependence on the 
feedback phase. The amplitude of the frequency shift is 
half of that of the change in damping, consistent with the 
discussion above. The phase where the damping is max- 
imized, coincides with the lowest resonator temperature 
and zero frequency shift. At this phase the system delay 
is compensated and a pure velocity-proportional feedback 
is applied to the resonator (i.e. ZAGfb — —n/2). The 
phase dependencies can also be calculated without any 
free parameters by using the values from the network 
characterization (Fig. |l]d). Figure 2 show that these are 
in good agreement with the feedback results. 
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To further cool the resonator, the feedback gain is in- 
creased at the optimal phase as indicated in Fig. ||. First 
the resonator temperature decreases rapidly with increas- 
ing gain due to the increased damping rate. However, by 
increasing the gain further, more of the imprecision noise 
u n is fed back as force noise. This causes a steady in- 
crease in for large gfb- The minimum temperature 
that can be reached is set by S UnUn and the solid line 
shows the predicted curve for velocity-proportional feed- 
back 0, |l3| calculated with the experimental parame- 
ters. The achieved minimum of 14.3 mK is close to the 
predicted lowest temperature of 14.0 mK. Note, that a 
temperature of 14.3 mK corresponds to an average ther- 
mal phonon occupation of n ~ ksTj^/hfu = 138 for a 
2 MHz resonator. The heterodyne DSP-based technique 
employed in this work thus successfully reaches the low- 
est temperature possible for the standard fully-analog ap- 
proach, but now applied to a high-frequency resonator. 

Another advantage of DSP-based feedback is that the 
transformation of v s and v c to a c and a s can be designed 
with almost arbitrary transfer characteristics, allowing 
implementation of optimal control strategies Jl4| . In the 
measurements in Fig. |^, the input signal is digitally fil- 
tered using a Fourier transform filter which is centered 
around fo and a tunable filter bandwidth Af. The fil- 
ter reduces the bandwidth of the feedback which pre- 
vents excess signal output outside the resonator band- 
width that may overload the detector or the amplifiers. 
The inset of Fig. || shows the root-mean-square output 
voltage as a function of gfb for three values of Af. For 
the full bandwidth (f s /2 — 410 kHz) an instability oc- 
curs around gfb — 0.23, which affects the cooling. The 
10 kHz bandwidth, which is used for the cooling curve 
of Fig. ||, has two orders of magnitude less actuation 
compared to the full bandwidth, enabling efficient cool- 
ing without affecting the closed-loop response as long as 
Af ^> 7ft/27r. This again illustrates the versatility of our 
digital quadrature feedback cooling platform. 



FIG. 3: Resonator temperature as a function of feedback gain, 
showing both the feedback measurements (symbols) and cal- 
culations for pure velocity-proportional feedback (solid line). 
The inset shows the root-mean-squared value of the actuation 
signal (a/ (of) + (of)) as a function of gain for different filter 
band widths. 
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